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Abstract 

We discuss renormalization of the non-relativistic three-body problem with 
short-range forces. The problem is non-perturbative at momenta of the order 
of the inverse of the two-body scattering length. An infinite number of graphs 
must be summed, which leads to a cutoff dependence that does not appear 
in any order in perturbation theory. We argue that this cutoff dependence 
can be absorbed in one local three-body force counterterm and compute the 
running of the three-body force with the cutoff. This allows a calculation of 
the scattering of a particle and the two-particle bound state if the correspond- 
ing scattering length is used as input. We also obtain a model- independent 
relation between binding energy of a shallow three-body bound state and this 
scattering length. We comment on the power counting that organizes higher- 
order corrections and on relevance of this result for the effective field theory 
program in nuclear and molecular physics. 



I. INTRODUCTION 



The three-body system provides non-trivial testground for ideas developed in two-body 
dynamics, and as such has a long and venerable history It is in general of very difficult 
solution, but when all particles have momenta much smaller than the inverse of the range of 
interactions it simplifies considerably while still retaining some very rich universal features 
[.y-@|. At such low energies, the two-body system can be attacked with many different 
techniques: effective range expansion, boundary conditions at the origin, etc.; a particularly 
convenient method for the extension to a many-body context is that of effective field theory 
(EFT) [||. Here we use EFT to solve the three-body system with short-range interactions 
in a systematic momentum expansion. 

Generically the sizes of bound states made out of A particles of mass m are all comparable 
to the range R of interactions. Similarly, the dimensionful scattering parameters — the two- 
body scattering length a 2 , the two-body effective range parameter r 2 , and so on — are also 
of the same order, R ~ a 2 ~ r 2 ~ . . . The A-body amplitudes reduce at low energy k 2 /m 
to perturbative expansions in kR ||. Less trivial is the case where the interactions are 
fine-tuned so that the two-body system has a shallow (real or virtual) bound state — that 
is, a bound state of size ~ a 2 much larger than the range of the interactions R ~ r 2 . This is 
a case of interest in nuclear physics, where the deuteron is large compared to the Compton 
wavelength of the pion l/m n , and in molecular physics, where shallow molecules such as 
the 4 He "dimer" can be over an order of magnitude larger than the range of the interatomic 
potential. 

In this fine-tuning scenario, the A-body system still can be treated in a perturbative 
expansion in fca 2 in the scattering region where all momenta are of order k <C l/a 2 . Bound 
states, however, correspond to k > a 2 and demand that a certain class of interactions 
be iterated an infinite number of times. In the two-body system, it can be shown ||[7| 
that one needs to sum exactly only two-body contact interactions which are momentum 
independent. This resummation generates a new expansion in powers of kR where the 
full dependence in fca 2 is kept. Contact interactions with increasing number of derivatives 
can be treated as perturbative insertions of increasing order. First and second corrections, 
for example, are determined by the insertion of a two-derivative contact interaction, which 
encodes information about the effective range; to these first three orders, all there is is 5-wave 
scattering with an amplitude equivalent (up to higher-order terms) to the first two terms 
in an effective range expansion. The EFT for the two-particle systems is thus equivalent to 
effective range theory it is valid even at k ~ l/a 2 and, in particular, for bound states 

of size ~ a 2 . 

There has been great progress recently in dealing with this problem in the two-body 
case 0. Ultraviolet divergences appear in graphs with leading-order interactions and their 
resummation contains arbitrarily high powers of the cutoff. A crucial point is that this cutoff 
dependence can be absorbed in the coefficients of the leading-order interactions themselves. 
All our ignorance about the influence of short- distance physics on low-energy phenomena is 
then embodied in these few coefficients, and EFT has predictive power. 

The question we want to answer in this paper is whether the leading two-body inter- 
actions are sufficient to approximately describe the three-body system in the same energy 
range, or whether we also need to include three-body interactions in leading order. The 
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answer hinges on the relative size of three-body interactions, so this question is intimately 
related to the renormalization group flow of three-body interactions with the mass scale 
introduced in the regularization procedure. This flow, in turn, depends on the behavior of 
the sum of two-body contact contributions to three-body amplitude as function of the renor- 
malization scale, or equivalently, as function of the ultraviolet cutoff A. What makes this 
problem different from standard field theory examples is that it does not have a perturbative 
expansion in a small coupling constant, and thus from the start involves an infinite number 
of diagrams. 

The extension of the EFT program to three-particle systems in fact presents us with a 
puzzle [jnj. The two-nucleon (NN) system has a shallow real bound state (the deuteron) in 
the 3 Si channel (and a shallower virtual bound state in the 1 S'o channel). Information about 
the three-nucleon system is accessible in nucleon-deuteron (Nd) scattering, which proceeds 
via two channels of total spin J = 3/2 and J = 1/2. Because of the Pauli principle we 
expect smaller three-body contact forces in the quartet channel than in the doublet channel. 
Assuming three-body forces to have a size completely determined by R according to naive 
dimensional analysis, the J = 3/2 Nd amplitude only receives contributions of three-body 
forces at relative 0((kR) 6 ), so up to (and including) relative 0((kR) 2 ) the J = 3/2 Nd 
amplitude depends on the AW amplitude only OJT2]. At low energies, the J = 3/2 Nd 
amplitude can be obtained by solving a single integral equation. We have shown that it 
is ultraviolet convergent, and its low-energy end independent of the cutoff; moreover, with 
parameters entirely determined by AW data, it predicts low-energy phase shifts in excellent 
agreement with the experimental scattering length |Tl| and with an existing phase-shift 
analysis fl2 |. The J = 1/2 Nd amplitude constructed out of the AW amplitude to the 
same order can be obtained from a pair of coupled integral equations. However, despite 
describing a sum of ultraviolet finite diagrams, numerical experimentation shows that it 
does not converge as the cutoff is increased. A system of three bosons exhibits a similar 
problem in the easier context of a single integral equation. For simplicity, we limit ourselves 
here to the latter case, which is of relevance to 4 He molecular systems. Similar arguments 
but more cumbersome formulas apply to the J = 1/2 three-fermion amplitude, which we 
will address in a later publication. 

Our study concerns momenta such that all forces can be considered short-ranged, but it 
is also relevant to higher momenta where we start to resolve the short-range dynamics. In 
general, what we have been considering as short-range dynamics will itself have structure 
and consist of longer- and shorter-range components. The latter will still be described by 
contact interactions in the new EFT appropriated to these higher energies. Elements of 
the discussion presented here will permeate this more complicated renormalization scenario. 
In the nuclear case, as we increase energy we start seeing effects of pion propagation, and 
interactions have both pion-exchange and contact components of similar sizes JL3]. It has 
recently been argued J7J that, because of fine-tuning, at moderate energies contact interac- 
tions are actually larger than pion effects by a factor of 2 or 3, so that pion interactions 
can be treated as corrections. A number of applications show that deuteron physics can be 



fruitfully dealt with this way |L4|. The momentum where this power counting breaks down 



is currently matter of controversy Our results hold unchanged throughout the region 
of validity of this power counting. 

After setting up our formalism in Sect. |TJ, we show in Sect. |TTT| that in general the sum 
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of two-body contact contributions to the particle/bound-state scattering amplitude, while 
finite, does not converge as the cutoff A is increased. We then present evidence in Sect. [IV 
that the leading-order cutoff dependence can be absorbed in a local three-body force with 
one single parameter A*. In particular, we derive an approximate analytic formula for the 
dependence of the bare three-body on the cutoff A. As a consequence, although input from 
a three-body datum is necessary, the EFT retains its predictive power. For example, if the 
three-body scattering length a 3 is used to determine A*, then the amplitude at (small) non- 
zero energies can be predicted and is cutoff independent. If A* is such that there exists a 
(ground or excited) bound state large enough to be within the range of the EFT expansion, 
its binding energy B 3 can be predicted as well. As an example, we consider the 4 He trimer. 
In any such physical system, A* is determined by the dynamics of the underlying theory 
as a certain (possibly complicated) function of its parameters. Different models for the 
underlying dynamics that are fit to the same two-body scattering data will correspond to 
different values for CI3 and B 3 , and in principle could cover the whole plane B 3 x 03. The 
EFT, however, predicts that in leading order these models can be distinguished by only 
one number, A*, and therefore that they would lie on a curve B 3 = £3(03). Indeed, in a 
nuclear context this has been observed, and is called the Phillips line 0. We discuss bound 
states and construct the corresponding bosonic line in Sect. M. In Sect. [VT| we discuss how 
corrections to this leading order can be handled in a systematic fashion. In Sect. |V11| we 
offer conclusions on our extension of the EFT program to three-particle systems with large 
two-body scattering lengths. Some of the results presented here were briefly reported in Ref. 

m 



II. LAGRANGIAN AND INTEGRAL EQUATION 

Particles with momenta much smaller than their mass m propagate as non-relativistic 
particles. If momenta are also small compared to the range of interaction R, the (bare) 
interaction can be approximated by a sequence of contact interactions, with an increasing 
number of derivatives. This is true regardless of the fine details of the interaction; infor- 
mation about these details is encoded in the actual values of the coefficients of the contact 
interactions. The most general Lagrangian involving a non-relativistic boson if) and invariant 
under small- velocity Lorentz, parity, and time-reversal transformations is 

c = ^(id + ^)ip - - ^(vW + . . . , (1) 

where the ellipsis stand for terms with more derivatives and/or fields; those with more 
fields will not contribute to three-body amplitudes, while those with more derivatives are 
suppressed at low momenta. 

The scope of this EFT in the two-body sector is by now well understood Q. The 
two-body amplitude will contain, in addition to "analytic" terms from two-particle contact 
interactions, also their iteration, which produces loops whose non-analytic terms are respon- 
sible for the unitarity cut. In leading order, one has to consider only the Cq term iterated to 
all orders, which is equivalent to the effective range expansion truncated at the level of the 
scattering length. The ultraviolet divergences found can all be absorbed in the parameter 
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Cq. The renormalized value of this parameter is Co = AiTa2/m. For 02 > there is a real 
S- wave bound state at energy — B 2 = —l/ma\. It can be further shown that the renormal- 
ized on-shell two-body amplitude forms an expansion in momenta that breaks down only 
at momenta of 0(1/ R) 0. In the fine-tuning scenario we are interested in, for momenta 
of magnitude typical of the shallow S*-wave bound state, k ~ l/a>2, the expansion is in the 
small parameter R/a 2 <C 1. We want here to extend this approach to the three-body system 
and establish the order in the EFT expansion in which a three-body force first appears. 



It is convenient |T7| to rewrite this theory by introducing a dummy field T with quantum 
numbers of two bosons, referred to from here on as "dimeron" (in analogy to dibaryon in the 
nuclear case). It is straightforward to show — for example, by a Gaussian path integration — 
that the Lagrangian Eq. ([!]) is equivalent to 

V 2 n 

£ = ^(id + + AT f T - ^=(7^ + h - c -) + hT^T^ip + ... (2) 

The arbitrary scale A is included in (H) only to give the field T the usual mass dimension 
of a heavy field. Observables will depend on the parameters appearing explicitly in Eq. @ 
only through the combinations g 2 / A = Cq and —3hg 2 /A 2 = D . 

Two-particle contact interactions are thus replaced by the s-channel propagation of the 
dimeron field; the two-body amplitude is completely determined by the full dimeron self- 
energy, that is, by the dimeron propagator dressed by particle loops. Three-particle contact 
interactions are likewise rewritten as particle-dimeron contact interactions, and so on. 

Let us consider elastic particle/bound-state scattering. (Three-body bound states man- 
ifest themselves as poles at negative energy. Both the inelastic channel and three-nucleon 
scattering involve the same type of diagrams and can be obtained by a straightforward ex- 
tension of the following arguments.) We restrict ourselves to external momenta Q ~ l/a 2 
and energies Q 2 /m ~ l/ma 2 . In leading order in a low- momentum expansion, the (bare) 
dimeron propagator is simply a constant i/A, while the propagator for a particle of four- 
momentum p reduces to the usual non-relativistic propagator 

iS(p) = — - 2 — , (3) 

P° ~ L + ie 

with S = 0(m/Q 2 ). 

A few of the first (connected) diagrams in perturbation theory in the coupling constant 
g are illustrated in Fig. [l|. In second order in the coupling, there is a single, tree diagram, of 
0(mg 2 /Q 2 ). The second diagram is of fourth order and has one loop; it is power-counting 
convergent and of 0((mg 2 ) 2 / 'AttQA). The third diagram is of sixth order and has two- 
loops; it is logarithmically divergent in the ultraviolet cutoff A and of 0((mg 2 ) 3 ln(A/A)/ 
(47rA) 2 ). It is clear that increasing the order brings factors of mg 2 Q/A7cA, so the divergences 
get progressively worse. At this point, one might jump to the conclusion that removal of 
divergences would require an increasing number of three-body force counterterms, starting 
with h ~ ((m^ 2 ) 3 /(4vrA) 2 ) ln(A/A). 

However, there are other diagrams that appear in the same orders: besides the bare 
"pinball" diagrams of Fig. [I], we a ls° have those which contain insertions of nucleon loops 
in dimeron propagators, the first few of which are shown in Fig. All but one of the 
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FIG. 1. Bare pinball diagrams up to 0(<? 6 ). A single (double) line represents a particle (bare 
dimeron) propagator. 

diagrams of fourth and sixth orders, contribute to (off-shell) wave-function renormalization. 
The remaining sixth-order diagram is of the form of the fourth-order bare pinball. The extra 
loop contributes a linearly divergent integral. The divergent piece oc A can be absorbed in 
g 2 /A, and for cutoffs A ~ 1/R the convergent terms that go oc Q 2 /A and higher inverse 
powers of A are of the size of other already disregarded higher-order terms. This is just 
two-body renormalization at work in a three-body context: the net effect of a nucleon loop 
in the dimeron propagator is the remaining, non-analytic term of 0(mg 2 Q/ATrA). 

If we were interested in momenta much smaller than the typical bound-state momentum, 
this would be essentially the end of the story. Things are more interesting for Q ~ l/a2, 
in which case the coupling constant expansion is an expansion in mg 2 Q/AirA ~ 1. In the 
two-body subsystem, all graphs built up from leading interactions are equally important. 
Whenever a bare dimeron propagator appears, so do all bubble insertions in it. It is not 
difficult to show that the two-body problem can still be renormalized by absorbing in 
g 2 /A the factors of A from more insertions of nucleon loops. (From here on g 2 / A refers 
to the sum of bare and loop contributions.) But what happens with the divergences in the 
graphs of Fig. [T|? 

The three-body amplitude can now be reexpressed in terms of the graphs in Fig. [I] but 
with the dressed dimeron propagator shown in Fig. §, 

-A + ^^-mpo + ^-ze + ze 

substituted for the bare propagator. The extra ^/ in the denominator improves ultraviolet 
convergence: now pinball loops carry factors of mg 2 Q/A7c(A + Q). All diagrams go as 






FIG. 2. Pinball diagrams with bubbles up to 0(<? 6 ). 
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FIG. 3. Dressed dimeron propagator. 



1/A 2 for Q ~ A, i.e., are power counting finite. The resummation into a dressed propagator 
accounts for cancellations among divergences found in diagrams with bare propagators. One 
could expect, based on perturbation theory experience, that for A ^> l/a 2 the low-energy 
end of this three-body amplitude is to a good approximation cutoff independent, and that 
the three-body amplitude converges as A — > oo. 

Unfortunately, things are not so straightforward. Since all diagrams are of the same 
order, 0(mg 2 Q/47rA) ~ 1, the three-body amplitude is actually the solution of an integral 
equation. This equation, including the three-body force, is depicted in Fig. || We choose the 
following kinematics: the incoming particle and bound state are on-shell with four-momenta 
(k 2 /2m, —k) and {k 2 /4m — B 2l k), respectively. The outgoing particle and bound state are 
off-shell with four-momenta (k 2 /2m — e, —p) and (k 2 /4m — B 2 + e, p); the on-shell point has 
e = k 2 /2m — p 2 /2m and p = k. The total energy is E = 3/c 2 /4m — B 2 . Denoting the blob 
in Fig. [|with this kinematics by it(k,p,e), we have 

it(k, p, e) = -2g 2 iS(-k 2 /4m - B 2 + e,p+ k) + ih 

+A Jj^jZ iS ( k2 / 2m ~ £ ~ Qo, [-2g 2 iS(-k 2 /4m - B 2 + 2e + q ,p + q) + ih 

iA(k 2 /4m - B 2 + e + q ,q) it(k, q,e + q ). (5) 

Here A = 1 for the bosonic case we are considering. The same equation is valid in the case of 
fermions, with different values of A. In particular, for three nucleons in a spin J = 3/2 state 
under two-body interactions only, this equation holds with A = — 1/2 and h = 0. (For three 
nucleons in a spin J = 1/2 state a pair of coupled integral equations results, with properties 
similar to the bosonic case.) 

After performing the go integration we can set e = k 2 /2m—p 2 /2m and, defining t(k,p) = 
t(k,p, k 2 /2m — p 2 /2m), we have 

2mg 1 



hi 
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+ p 2 + mB 2 + p- k 



+8ttA 



d 3 q t(k, q) 



(2tt)3 _i_ + /3£ _ mE 

(12 V 4 



3k 2 
4 



1 h 



+ mB 2 + q 2 +p 2 +p ■ q 2mg 2 



(6) 



Projection on the S'-wave is obtained by integration over the angle between p and k, resulting 



in 



t(k 1 (p 2jr Vk + k 2 — mE\^^ 

pk \p 2 — pk + k 2 — mE I 
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T± ■ 
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FIG. 4. The amplitude T for particle/bound-state scattering as a sum of dressed pinball and 
three-body- force diagrams (first and second lines) and as an integral equation (third line). 



2A M 
+ — / dq- 



t{k,q)q 2 



n Jo _± + m _ mE 

a,2 V 4 

The scattering amplitude T{k) is given by 

T(k) = \fZt{k,k)yfZ, 

where Z is given by 



1 / p 2 + pq + q 2 — mE \ h 
In 



Z - 1 = i |_ (iA (P))- 



2 2 

m g 



P0 = ~B2 



It is customary to define the function a(k,p), 

a(k,p) mo 2 Zt(k,p) 



V 2 - k 2 Svr _± + J^-mE 

a.2 V 4 



that on shell reduces to 



m 



a(k,k) = —T(k). 



pq \p 2 — pq + q 2 — mE J mg 2 



In particular, a(0, 0) = —03. The equation satisfied by a{k,p) is 

2A f A q 2 

a(k,p) = M(k,p;k) + — / dq M(q,p;k)— -a(k,q), 

Ti Jo q z — k A — le 



(7) 



(8) 



(9) 



(10) 



(11) 



(12) 



with the kernel 
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M(q,p; k) 




1, / q 2 + qp + p 2 
— In 
pq 



mE 



q 2 — qp + p 2 — mE 



+ 



mg z 



(13) 



Eqs. ( |T2] ) and ( |I~3"D reduce to the expressions found in Refs. p^8|JT0H l2| when h = 0. Note 
that the perturbative series shown in Fig. |] corresponds to a perturbative solution of the 
integral equation for small A. 

The solution of Eq. (|12|) is complex even below the threshold for the breakup of the two- 
particle bound state due to the ie prescription. To facilitate our discussion we will use below 
the function K(k,p) that satisfies the same Eq. (|12"D as a(k,p) but with the ie substituted 
by a principal value prescription. K(k,p) is real below the breakup threshold. a(k, k) and, 
consequently, the scattering matrix can be obtained from K(k,p) through 



a(k, k) 



K(k,k) 
-ikK(k,k)' 



(14) 



III. THE PROBLEM 

In order to understand the ultraviolet behavior of the theory, let us first take h = 0. 
When p ^> l/a 2 (but k ~ 1 / a 2 ) , the inhomogeneous term is small (0(l/pa 2 )), the main 
contribution to the integral comes from the region q ~ p, and the amplitude satisfies the 
approximate equation 

K(k, P ) = « / A * In ( 4±H±4) K(k, q) . (15) 
v37r Jo q \q 2 -pq + p 2 J 

Now, in the limit A — ► oo there is no scale left. Scale invariance suggests solutions of the 
form K(k,p) = p s , which exist only if s satisfies 

8A sin^f 

1 - -= §- = 0. 16 

V3s cos f K J 

If K(k,p) is a solution, K{k,p 2 /p) is also a solution for arbitrary p+. Because of this 
additional symmetry, the solutions of Eq. ( |I~6"D come in pairs. 

For A < A c = ~ 0.4135, Eq. (O) has only real roots. For example, if A = —1/2, there 



are solutions s = ±2, ±2.17, ... In this case an acceptable solution of Eq. (|15 ) decreases 
in the ultraviolet, K(k,p ^> 1/0-2) = Cp~^. For finite cutoff, the solution should still have 
this form as long as p < A. The overall constant C = C(A) cannot be fixed from the 
homogeneous asymptotic equation, but is determined by matching the asymptotic solution 
onto the solution at p ~ l/a 2 of the full Eq. (0). Because of the fast ultraviolet convergence, 
the full solution to Eq. (|T^) is expected to be insensitive to the choice of regulator, so that 
a well-defined liniA^oo C(A) can be found. This behavior can indeed be seen in a numerical 
solution of Eq. (0) 

For A = 1, on the other hand, there are purely imaginary solutions: s = ±^So, where 
s ~ 1.0064. The solution of Eq. ([15]) as A -> 00 is 
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FIG. 5. Amplitude K(0,p) as function of the momentum p. Full, dashed, and dash-dotted 
curves are for H = and A = 1.0, 2.0, 3.0 x W^a^ , respectively. Dotted, short-dash-dotted, and 
short-dashed curves are for A = 1.0 x 10 4 a2 1 and H = —6.0, —2.5, —1.8, respectively. 



K(k,p ^> l/a 2 ) = C cos [ s ln— ] . 

V P*J 



(17) 



Again, for finite cutoff this should hold at least for p C A, with p* = p*(A) determined by 
the cutoff. On dimensional grounds, 



p*(A) = exp(-5/s ) A, 



(18) 



where 5 is a dimensionless, cutoff independent number. 

There are now two constants C(A) and 5: the solution of Eq. (|T2| ) is not unique in 
the limit A — > oo |19|]. The undetermined phase p* arises from the symmetry K(k,p) — > 
K(k,pl/p). When A = — 1/2 we can dismiss the solution that grows in the ultraviolet since 
the integral in the Eq. (|12"D will not converge, but in the A = 1 case there is no way to 
select a preferred oscillatory solution. The same problem exists for any A that yields purely 
imaginary solutions for s in Eq. fllTf ). From now on, for definiteness we specialize to the 
case of purely imaginary s. 

Solutions K(0,p) of Eq. flT2|) for A = 1 and finite A (but h = 0) are shown in Fig. [5[ We 
have verified that in the region l/a 2 p A the solutions indeed have the form (|17D , with 
So = 1.02 ± 0.01. The phase p± was found approximately linear in the cutoff in accordance 
with Eq. (18|): 5 = 0.76±0.01 is cutoff independent. We have also found that the amplitude 



is 



C 



7 



cos (s ln(p*a 2 ) + e) 



(19) 



with 7 = 1.50 ± 0.03 and e ~ 0.08, both cutoff independent. This generates the striking 
meeting points at p n = a^expd^n + l)ir — e]/s ) ~ 0.92a2 1 (22.7)^ n+1 - ) , n an integer, that 
can be seen in Fig. [5| 
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FIG. 6. Amplitude K(0,p) for ^,1,2,3, in the family of Ao = 5.7202 1 , corresponding to 
H(A n ) = 0, for a 3 = -2.0a 2 . 

Apart from these meeting points, the amplitude does not show signs of converging as 
A — ► 00. We are forced to conclude that terms that are 0(Q/A) cannot be dropped, as it is 
usually done in calculations in finite orders in perturbation theory. The presence of these 
terms determines a unique phase in the asymptotic region. The solution for small p is to 
be matched at an intermediate scale to the large-p solution, so the cutoff dependence leaks 
into the small-p region. Small differences in the asymptotic phase lead to large differences 
in, for example, the particle/bound-state scattering length. 

Since K(k,p ^> 1/02) is given by Eqs. (fL7l) -([L9l) , it is the same within a discrete family 
of cutoffs A n = A exp(wr/so) — A (22.7) n , n an integer. We expect that the low-energy 
solution will then be invariant for cutoffs in this family. If A is fixed in such a way that 
the three-body scattering length a% is reproduced, then the same a 3 should result from any 
of the A ra 's in the same family. Moreover, the amplitude K(k,p) around p = k = should 
be similar as well. This is indeed what we find: in Fig. || we plot the low-energy solutions 
for cutoffs in the family of A = 5.72a2 1 (Ai = 122. 2a2 1 — 22.7A , etc.) fitted to give 
a 3 = — 2.0a 2 (cf. Fig. g). For the appropriate a 3 /a 2 , we find the same cutoffs as in Ref. [ED| . 



This generalizes the results of Ref. [pOf] : we see that a family of cutoffs determines not only 
the scattering length but also the (low-)energy dependence of the amplitude. As we vary 
families (i.e., Ao), however, the low-energy behavior of the amplitude shows strong cutoff 
dependence. This leakage of high-momentum behavior into the low-momentum physics is 
indication that we are not performing renormalization consistently with our expansion. 

Note that if one were to truncate the series of diagrams in Fig. [| at some finite number 
of loops one would miss the asymptotic behavior of K(k,p) that generates this cutoff de- 
pendence. This is because sq (and its expansion in powers of A) vanish in a neighborhood 
of A = 0. The truncation of the series in Fig. |] is equivalent to perturbation theory in A, 
and cannot produce a non-vanishing sq. We are here facing truly non-perturbative aspects 
of renormalization. 
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IV. THE SOLUTION 



Faced with such a problem, the only way to eliminate this cutoff dependence is to modify 
our leading order calculation, that is, to change our accounting of the higher-energy behavior 
of the theory through the addition of at least one new interaction. 

In order to do so, we can follow essentially two routes. One is to revise our systematic 
expansion in the two-body system. The task in this case is to enlarge the EFT in order to 
incorporate higher-energy degrees of freedom that have so far been treated only in something 
akin to a multipole expansion. Generically, the next mass scale above l/a 2 is l/r 2 , and 
incorporating momenta Q ~ l/r 2 means that amplitudes certainly can no longer be treated 
in a Qr 2 expansion. With these new elements, a new small parameter has to be identified. 
The resulting, new leading-order two-body amplitude might very well exhibit sufficiently nice 
ultraviolet behavior to guarantee the disappearance of the pronounced cutoff dependence 
found above in the three-body amplitudes. 

While this is a perfectly legitimate way to proceed, it is in most cases not a simple task. 
It demands a qualitative new understanding of the shorter-range dynamics of the system 
under study. It is also specific to that system since the dynamics at Q ~ l/r 2 will have to 
be treated in more detail; its consequences will not be common to all systems with a large 
a 2 . Moreover, the modifications of the leading-order two-body amplitude, if any, might not 
resolve the leakage problem in the three-body system. We follow here the other possible 
route. By construction, this approach preserves the expansion of the two-body amplitude; 
as such it is simple, model independent, and designed to function. 

The observed cutoff dependence comes from the behavior of the amplitude in the ultravi- 
olet region, where the EFT Lagrangian, Eq. (Q), is not to be trusted. When the low-energy 
expansion is perturbative the cutoff-dependent contribution from high loop momenta can be 
expanded in powers of the low external momenta and cancelled by terms in the Lagrangian, 
since those also give contributions analytic in the external momenta. Thus all uncertainty 
coming from the high momentum behavior of the theory is parametrized by the coefficients 
of a few local counterterms. The present case is complicated by the fact that the cutoff 
dependence of the amplitude is non-analytic around p = 0. This, however, does not mean 
that the renormalization program in this low-energy EFT is doomed: a three-body force 
term of sufficient strength contributes not only at tree level, but also in loops dressed by 
any number of two-particle interactions. We want to show that the bare three-body force 
coefficient h(A) can be combined with the the problematic parameter p*(A) in order to pro- 
duce an amplitude K(k,p) which is cutoff independent, at least around p — k — 0. Because 
the cutoff dependence generated by the two-body amplitude is a somewhat complicated 
oscillation, we expect a similar behavior of the h(A) that ensures cutoff independence. In 
particular, h(A) must be such that it vanishes for "critical" cutoffs A n that belong to the 
family that gives the desired 03. For any other cutoff, h(A) must be non- vanishing. 

So, let us now consider the effects of a non- zero h(A). It is convenient to rewrite the 
three-body force as h(A) = 2mg 2 H (A) / A 2 . In the remainder of this paper we will refer 
to H(A) as the bare three-body force. H(A) has to be at least big enough so as to give a 
non-negligible contribution in the p ~ k ~ A region. This means that the dimensionless 
quantity H(A) has to be at least of 0(1). We will look for a solution initially assuming a 
"minimal" ansatz H(A) ~ 1; we subsequently show via an unbiased numerical analysis that 
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this assumption is justified. 

Such a three-body force has the feature that its contribution to the inhomogeneous term 
is small compared to the contribution from the two-body interaction, as it is at most p/A 
of the latter. When p ^> l/a 2 (but k ~ l/a 2 ) , the amplitude satisfies a new approximate 
equation 



K{k,p) 



4A r A dq 
Jo q 



In 



q 2 + pq + p 2 ' 
q 2 — pq + p 2 



+ 2H(A) 



pq_ 

A 2 



K{k,q). 



(20) 



For p ~ A the term proportional to H(A) becomes important and the solution K(k,p ~ A) 
has a complicated form. In the range l/<22 <C p <C A, on the other hand, the three-body force 
term is suppressed by p/A compared to the logarithm and can be disregarded. Consequently, 
the behavior ([17]) is still correct in the intermediate region. The effect of a finite value of 
H can be at most to change the values of the amplitude C and the phase 5, which become 
dependent on H. As shown in Fig. |5], this feature is confirmed by numerical solutions: while 
different values of the three-body force H (at a fixed A) preserve the form of the solution, 
the phase (and amplitude) are changed. If H is chosen to be a function of A in such a 
way as to cancel the explicit A dependence, we can make the solution of Eq. fll2T) cutoff 
independent for all p C A. In particular, the scattering amplitude that is determined by 
the on-shell value K(k, k) with k ~ l/ci2 will be cutoff independent as well. For this to be 
possible C and S must depend on the same combination of A and H. Thus H(A) must be 
chosen such that 



- s In A + 5(H(A)) = -s hiA* 



(21) 



where A* is a (cutoff independent) parameter fixed either by experiment or by matching with 
a microscopic model. Eq. ( pT|) simply means that p*(/i(A),A) = A* is cutoff independent, 
and the l/a2 CpC A solution is 



i^(/c,A>p> l/a 2 ) 



7 



V 



cos(s m(A*a 2 ) + e) V A 



(22) 



Matching with the p ~ l/a2 solution should then determine the scattering length a 3 = 03 (A*) 
and the low-energy dependence of the amplitude. 
For p ~ 1/02, Eq. ( jl2]) becomes 



K(k,p) 



0,2 V 4 



1 fp 2 +pk + k 2 - ME\ 
pk n Vp 2 -pk + k 2 - ME ) 

2\ ( p 2 +pq + q 2 - ME\ qK(k,q) 
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+ — / dq 
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pq + q 2 — ME J q 2 — k 2 



q- 



H(A) 
A 2 



K(k,q), 



(23) 



where /x is an arbitrary scale such that /i< A, and we have dropped terms which are smaller 
by powers of p/A or /i/A. We want K(k,p ~ l/a 2 ) cutoff independent and determined 
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(a) 



(b) 



FIG. 7. Inhomogeneous terms of the integral equation 
kernels. 



2|): (a) two-body and (b) three-body 



completely by A*. Then all cutoff dependence contained in the //-to- A integral — in the 
integration limit, implicit in K(k,p ^> 1 / a 2 ), and in H(A) — has to combine to produce a 
A-independent result. Achieving this means that K(k,p ~ 1/02) will depend only on A* not 
just at p = but also for a range of p's, because the high-momentum contributions from 
particle exchange and from the three-body force have the same momentum dependence. The 
high-momentum modes in Fig. 0(a) can be absorbed in Fig. 0(b). 

We can now obtain an approximate expression for H(A). For this purpose, let us consider 
Eq. for two different values of the cutoff A and A' > A, whose solutions we denote by 
K(k,p) and K'(k,p), respectively. The equation for K'(k,p) will have the same form as the 
one for K(k,p) except for some extra terms: 




- + 



g(A) 
A 2 



K'(k,q) 



' H(A') H{A)\ 



A' 2 



A 2 



dq K'(k, q) 



Assuming that K'(k,p) has the same phase cos(so ln(p/A*)) as K(k,p) even for p 
(E24T) becomes 



(24) 



A', Eq. 



4A 1 

{ — [sin(s ln(A7A*) - arctg(l/s )) + H{A') sin(s ln(A'/A*) + arctg(l/s ))] 



nJl + s 2 " A 



"(A' -A)}, 

plus terms terms that are smaller by p/A. We can make the terms in fl25|) vanish if 



H(A) 



sin(s ln(£) - arctg(^)) 



sin(s ln(£) + arctg(^)) 



(25) 



(26) 



Note that H(A) is periodic in A: H(A n ) = if (A) for A n = Aexp(n7r/s ) ~ A(22.7) n , n 
an integer. In particular, the three-body force (|26|) vanishes at the critical cutoffs A n = 
A exp(nn / s ) with A = A^ exp(arctan(l/s )), as anticipated. Clearly, H is periodic in A* 
as well. 

Since with such H(A) the inhomogeneous equation for K ; (k,p) is nearly the same as the 
equation for K(k,p), it follows that K'(k,p) has also the same amplitude as K(k,p) in the 
intermediate region. That is, H(A) chosen like Eq. (|26| ) exactly compensates any change in 
cutoff, so that K'(k,p) = K(k,p) for all values p< A (up to terms suppressed by p/A). As 
a consequence, the on-shell A'-matrix K(k, k) will be A independent as long as k A. 
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FIG. 8. Three-body force H as a function of the cutoff A for a 3 = 1.56a2: numerical solution 



(dots) and Eq.(§D witn A * = 19.5o2 (solid line). 



In order to check these arguments we have numerically solved Eq. (|12|) for the amplitude 
K(k,p) with a nonvanishing H(A), in the case A = 1. We have fixed the scattering length 
a 3 = —K(0, 0), and for several cutoffs we have determined the three-body force that is neces- 
sary to keep a 3 unchanged. With such a procedure, we have obtained H(A) for different A's. 
In Fig. |8| we plot the corresponding numerical values for the case a 3 = 1.56a 2 , together with 
our approximate analytical expression (|26|) with A* = 19. ^a^ 1 . The remarkable agreement 
confirms that the solution for H(A) is very close to the minimal one, all the way up to very 
high values of A. Note that the three-body force vanishes at the critical cutoffs in the family 
generated by A = A* exp(arctan(l/s )) = 42.6a2 1 .[] 

In Fig. |9| we show results for the corresponding K(k, k)~ x = kcotS, where S is the 
5"- wave phase shift for particle/bound-state scattering, for different cutoffs in the case a 3 = 
1.56a2- Comparison with Fig. |5] shows that by introduction of the three-body force we have 
succeeded in removing the cutoff dependence of the low-fc region, thus generalizing Fig. ||| 
to cutoffs other than critical. As argued above, k cot 5 is insensitive to A as long as k <C A. 
The effective range, for example, is predicted as r 3 = 0.57ci2. In Fig. [X0| we plot the 5-wave 
phase shifts directly. 

Similar results can be obtained for other values of a 3 . By repeating this procedure we 
can in fact obtain a 3 = a 3 (A*) or vice-versa. This is shown in Fig. |ll|. The periodicity 
of the three-body force with A* is apparent. If the underlying theory is known and can be 
solved, then A* can be determined in terms of underlying parameters, and a 3 can also be 
predicted. Otherwise, we can use a 3 or any other low-energy three-body datum to determine 
A* empirically. In either case, once the two-body scattering length 02 is fixed, the low-energy 
particle/bound-state scattering amplitude is completely determined by the mass scale A* 
contained in the three-body force h. 



Note that the critical cutoffs in Fig. ^ were obtained for a different a 3 /a2 . 
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FIG. 9. Energy dependence for A* 
(A = 42.6, 100.0, 230.0, 959.0 x a 2 l ). 
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k cot 5 as function of k for different cutoffs 



V. BOUND STATE 



We can extend the preceding analysis to the three-body bound-state problem. In this 
case the relevant equation to be solved is Eq. (|l2|) without inhomogeneous terms. The latter 
did not play an important role in our ultraviolet arguments, so the same arguments hold for 
any bound state with binding energy B 3 comparable with 1/maf, i.e., for any bound state 
with the dimensionless binding energy 63 = ma\B^ = 0(1). Such bound states are shallow, 
having a size ~ \j\pmB~l = a<il\fb% comparable to a 2 , and thus should be within range of 
the EFT. In principle, all bound states with size larger than ~ r 2 should be amenable to 
this EFT description. Their properties will be determined in first approximation by only 
a 2 and 0,3, or equivalently Co and A*, while more precise information can be obtained in an 
expansion in \fbj,r 2 ja 2 . 

In Fig. O we plot binding energies for a range of cutoffs, with the three-body force 
adjusted to give a fixed scattering length a 3 = 1.56a 2 . (With an appropriate a 3 , our results 
reduce to those of Ref. [^0| when the cutoff is at one of the critical values A n .) As we can 
see, for this value of A* there exists a shallow bound state at 63 ~ 2 whose binding energy is 
independent of the cutoff. This bound state has a size ~ 0.7ct2 and can thus be studied within 
the EFT. As we increase the cutoff, deeper bound states appear at A n = A exp(n7r/so) with 
n an integer and A — 10, so that for A n _! < A < A n there are n + 1 bound states. The 
next-to-shallowest bound state has 63 ~ 70 and thus a size ~ 0.1a2. If the underlying theory 
is such that r 2 < a 2 /10, a cutoff A > l/r 2 allows us to deal with the next-to-shallowest 
bound state. The next deeper bound state has 63 ~ 4 x 10 4 and thus a size ~ 0.01a2, and 
so on. These most likely lie outside the EFT. 

Properties of the shallowest bound state, if within the EFT, are model independent and 
thus the most interesting. For a fixed a 2 , they are determined by A*. In Fig. [13] we plot 
the binding energy B3 as function of A*. This should be compared with the behavior of 
03 (A*) shown in Fig. [11]. It is clear that at A* ~ Aa 2 1 a bound state appears at zero 



03 



energy. As A* increases, this state gets progressively more bound until at A* ~ 22.7 x 4a 2 
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FIG. 10. S"-wave phase shifts 5 as function of energy E, as predicted by the EFT when 
a 3 = 1.56a 2 (A = 42.6, 100.0, 230.0, 959.0 x a 2 x ). 

a new shallowest bound state appears; the picture repeats indefinitely. Whenever a bound 
state is close to zero energy, a 3 is large in magnitude, negative if the bound state is virtual 



and positive if real. This can be seen in Fig. 11 



Eliminating A*, we obtain an universal curve B 3 = .63(03), as in Fig. |14[ Qualitative 



features of this curve can be understood from Figs. |TT| and |13[ For a 3 large in magnitude, 
small variations of A* lead to large variations of a 3 but small variations of B 3 , so the curve 
flattens out at both ends. For 03 large and positive, there is a shallow real bound state; for 
03 large and negative, there is a shallow virtual bound state but the curve tracks the deeper 
real bound state. For 03 ~ a 2 , a 3 and B3 vary with A* at a similar rate, and the curve 
interpolates between the two ends. 

This curve is known in the three-nucleon the Phillips line [B| . It has been derived 

before in the context of models for the two-particle potential that differed in their high- 
momentum behavior. Varying among two-particle potential models one could expect to fill 
up the B3 x a 3 plane. In the EFT, the high-momentum behavior of the two-particle potential 
is butchered, which causes no trouble in describing two-particle scattering at low-energies, 
but — in the peculiar way described here — does require a three-body force. It is crucial 
that the EFT does retain a systematic expansion, so that in leading order it requires only 
one local three-body counterterm, determined by A*. Varying among two-particle potential 
models is thus equivalent to varying the one parameter A* of the three-body force. This 
spans a single curve B 3 = -83(03) in the B 3 x a 3 plane. Our argument suggests that the 
Phillips line is a generic phenomenon and provides a simple explanation of its origin. 



VI. HIGHER ORDERS 



The corrections to the calculations in the previous sections come from operators not 
explicitly written in Eq. (jy) and are suppressed by powers of kR or R/a 2 . The first correction 
comes from terms involving two derivatives and four nucleons fields ||: 
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FIG. 11. Particle/bound-state scattering length 03 as function of the three-body force param- 
eter A*. 



Co 



— (^v^ f vV + •••)• 



(27) 



They account for the effective range term in the effective range expansion of particle-particle 
scattering, and one finds C 2 = 7ra 2 r 2 /m. As it was done before, ( pTj) can also be generated 
by integrating out the T field from Eq. (|2|) if the extra term 



(28) 



where = 4771C2/C0, is added to Eq. (0). Indeed, integration over the auxiliary field T 
generates, besides the terms shown in Eq. ([!]), also a four-nucleon term involving two space 
derivatives or one time derivative. As usual, a nucleon-field redefinition involving the equa- 
tions of motion can be performed to trade the time derivative by two space derivatives. This 
field redefinition does not change on-shell amplitudes. As far as observables are concerned, 
adding fl2"8|) to Eq. (0) is thus equivalent to adding (|27|) to Eq. ([!]). 

The correction proportional to r 2 is then given by the diagram in Fig. [K|(a) (plus 
corresponding wave function renormalization pieces) and includes one insertion of the kinetic 
term, Eq. (|2"8"|). The contribution of this graph to t(k,p) is 
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As discussed above, t(k,p) behaves for large p as t(k,p) ~ A(k) cos(s log(p/A*))/p, so the 
diagram is naively linearly divergent. The identity 
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FIG. 12. Three-body binding energies B% as functions of the cutoff A for A* = 19.5a 2 



can be used to rewrite Eq. (|29|) and it shows that the most divergent piece (the one due 
to the constant on the right hand side of Eq. (|30D) actually vanishes. This is because 
in this term the integral can be closed in the complex plane without circling any pole. 
The only contribution comes from the second term in Eq. P0| ) but this term is further 
suppressed in the ultraviolet, resulting in a logarithmically divergent contribution from the 
diagram in Fig. |T5|(a). This divergent piece has a complicated dependence on the external 
momentum k since it is proportional to A 2 (k) (the dependence on k of the other terms 
included in the graph is unimportant in the ultraviolet). This is an unusual situation. In 
perturbative calculations the divergent part is simply a polynomial in the external momenta 
and it can consequently be absorbed in a finite number of local counterterms appearing 
only at tree level. Here the dependence on the external momenta is more complicated but 
the counterterm appears in graphs with an arbitrary number of leading-order interactions. 
To see this, let us split the three-body force coefficient h into a leading order piece 
(the same considered in the previous sections) and a sub-leading piece that will be 
included perturbatively. The inclusion of the sub-leading three-body force proportional to 
foi 1 ) generates the three diagrams shown in Fig. |15|(b,c,d). Simple power counting shows 
that the graph in Fig. |l5|(d) is the most divergent of the three. The important point is 
that the external momentum dependence of the most divergent part of the graph in Fig. 
151(d) is given by A 2 (k), again because the dependence on the external momentum k of the 



nucleon propagators in Fig. |T5|(d) can be discarded in the ultraviolet region. Since the k 
dependence of the divergent part is the same in Fig. 0(a) and Fig. P(d), can be 
chosen as a function of the cutoff so as to make the sum of these two graphs finite. As 
always, the renormalization procedure leaves a finite piece in h undetermined that should 
be fixed through the use of one experimental input. This conclusion seems to agree with 



Ref. [21]. More details of the calculation of the range correction as applied to the case of 



three nucleons in the triton channel are left for a future publication. 
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FIG. 13. Binding energy B3 of the shallowest bound state as function of the parameter A*. 
VII. DISCUSSION AND CONCLUSION 

Our results should hold for 4 He atoms. In fact, it has recently been established that the 
two- 4 He bound state ("dimer") is very shallow, with an average size (r) = 62 ± 10 A m, 
more than an order of magnitude larger than the range of the interatomic potential. The low- 
energy two- 4 He system should then be describable by contact two-body forces. In leading 
order, the measured size translates into a scattering length a 2 = 124 A, which determines 
the strength of the contact interaction. Unfortunately, although the three- 4 He ("trimer") 



has been observed [p3l , there seems to be no low-energy information on its properties nor 
on 4 He/dimer scattering. At least one three-body datum is needed to determine A* and use 
the EFT to make predictions. 

Until such datum becomes available, we can only illustrate the method by using a phe- 
nomenological 4 He- 4 He potential as a model of a microscopic theory. We select a potential 



24[ which is consistent with the recent measurement of the dimer binding energy. It gives for 
the two-body system a 2 = 124.7 A and r 2 ~ 7.4 A. Three-body calculations are much more 
difficult to perform with such a phenomenological potential. An estimate for the 4 He/dimer 
scattering length is a 3 = 195 A. Ground and excited bound states have been reported; es- 
timates for the shallowest bound state place it in the range B3 = 1.04 — 1.7 mK, while a 
deeper state lies around B% = 0.082 — 0.1173 K. There exist a prediction for the low-energy 
S'-wave phase shifts, albeit for a different potential, but r% could not be determined. 

Using such model we can estimate the range of validity of the EFT in momentum to be 
~ l/r 2 ~ 0.14 A -1 , and the leading order to give an accuracy of ~ r 2 /a2 — 0.06, or about 
10%, at momenta Q ~ l/a2- Using this potential's 03/02 = 1.56 we can determine A*, and 
predict both the energy dependence in 4 He/dimer scattering and the trimer binding energy. 
In fact, in Figs. [], [II], and [EJ we have used exactly this value of 03. As a consequence, 
Fig. H represents our lowest-order prediction for kcotS for atom/dimer scattering, from 



which we can for the first time extract an effective range r% = 71 A. Fig. [10] displays the 



S-w&ve phase shifts themselves, and is in qualitative agreement with the result of Ref. |24 
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FIG. 14. Phillips line: binding energy B% of the shallowest bound state as function of the 
particle/bound-state scattering length 03. 

for a similar potential, being obtained here at a negligible fraction of computer time. From 
Fig. 12 we predict a bound state at B% = 1.2 111K, which is certainly within the EFT. The 
next-to-shallowest bound state is small enough to be at the border of EFT applicability. For 
a sufficiently large cutoff, we find B 3 = 0.057 K, but in the best case scenario corrections 
from higher orders should be a lot larger than 10%, and very important. 

Because of the similarity of the integral equations, our arguments should be relevant for 
systems of three fermions with internal quantum numbers as well. We are in the process of 
verifying this for the three-nucleon system in the J = 1/2 channel, where we will be able 
to check our predictions against the energy dependence of existing phase shift analyses of 
neutron-deuteron scattering and against the triton binding energy ||25|| . We also plan to 
extend our calculations to the four-and-more-body system and search for the existence of 
other leading few-body forces. 

These results will hold in an EFT where the pion has been integrated out, which is 



valid for Q 



< 



In an exciting new development, it has been argued in Ref. that 



there is a region of momenta above m n where, although pions have to be kept explicitly, 
their effects are sub-leading. The leading two-body operators are thus the same as in the 







(a) 



(b) 



(c) 



(d) 



FIG. 15. 0(r2/a2) diagrams: (a) one insertion of pq ) in the dimeron propagator, denoted by 
a cross; (b,c,d) one insertion of a correction to the three-body force h, denoted by a circle. 
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"pionless" theory. A number of examples seem to corroborate this picture JT3 . In this case, 



our leading-order results will be valid in this "pionful" theory as well, suggesting that our 
bound-state calculation will provide a reasonable estimate of the triton binding energy 

In conclusion, we have provided analytical and numerical evidence that renormalization 
of the three-body problem with short-range forces requires in general the presence in leading 
order of a one-parameter contact three-body force. This frames results obtained earlier with 
particular models within a larger, model-independent picture. It opens up the possibility 
of applying the EFT method to a large class of systems of three or more particles with 
short-range forces. 
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